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ABSTRACT 

A  new  technique  for  calculating  wave  functions  for  atoms 
and  ions  is  developed.   Formal  asymptotic  expansions  with 
largeness  parameter  Z   (nuclear  charge)  are  determined  through 
first  order  in   "^/Z-   In  so  doing  the  Schrodinger  Equation  for 
a  general  atomic  system  with  N  electrons  is  reduced  to  a 
set  of  first  order  partial  differential  equations  which  are 
solved  recursively.    Screening  and  correlation  are  exhibited 
explicitly  in  the  resulting  asymptotic  atomic  wave  functions. 
Applications  made  to  the  ground  state  of   2   electron  systems 
show  that  the  asymptotic  wave  function  obeys  the  vlrial  theorem 
through  first  order.   Magnetic  susceptibilities  within   5/0   of 

the  accepted  values  are  obtained  for  helium  and  singly  ionized 

2 
lithium.   Other  expectation  values,   l/2(s-|^  +S2  )?  s-^^        ^^'^ 

s-,p  are  found  to  be  in  excellent  agreement  with  a  variational 

calculation  utilizing  IO78  parameters.   In  determining  the 

asymptotic  wave  function  in  regions  about  certain  singular 

points  by  a  "matching  technique",  in   Z/z  terms  enter.   The 

omission  of  such  terms  in  ordinary  varlational-perturbation 

wave  functions  may  result  in  a  loss  of  accuracy  when  computing 

expectation  values  other  than  the  energy. 
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ASYMPTOTIC  WAVE  FUNCTIONS  FOR  ATOMIC  SYSTEMS 

1 .   Introduction. 

Recently  approximate  ground  state  solutions  to  the 
nonrelatlvlstlc  Schrodlnger  equation  for  many  electron  atoms 
have  received  considerable  attention.   Electron  wave  functions 
giving  a  reasonably  correct  description  of  the  electron  density 
throughout  configuration  space  would  facilitate  calculation 
of  weak  interaction  effects.   Of  late,  there  have  been 
refinements  in  the  experimental  measurements  of  magnetic 
susceptibility,  electric  polarlzability  and  quadrupole  coupling 
of  atoms.   The  utilization  of  more  accurate  wave  functions 
would  allow  calculation  of  these  effects  to  keep  pace  with 
experiment.   More  accurate  atomic  scattering  factors  could 
also  be  obtained.   In  addition,  if  accurate  excited  state 
wave  functions  were  known,  transition  probabilities  and  optical 
absorption  and  inelastic  scattering  cross-sections  could  be 
calculated  with  much  greater  precision. 

The  calculation  of  ground  state  atomic  wave  functions  is 
usually  performed  by  either  a  Hartree  (Hartree-Pock)  calculation 
or  by  some  other  variational  technique.   It  is  well  known  that 
Hartree  calculations  usually  give  the  energy  of  the  ground 
state  to  within  1  per  cent.    However  these  solutions  do  not 
incorporate  interelectronic  repulsion  correctly  since  they 
describe  each  electron  as  moving  in  the  average  central  field 
due  to  all  the  others  and  these  are  distributed  throughout 
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their  wave  functions  according  to  the  statistical  probability 
that  they  will  be  found  at  any  one  coordinate.   Hence  the 
electrons  do  not  avoid  one  another  because  of  their  mutual 
coulomb  repulsion  but  tend  to  avoid  regions  of  large  average 
potential  energy.   Also  such  techniques  are  derived  from  a 

variational  principle  and  their  application  to  excited  states 

2 
leads  to  problems  of  orthogonality. 

Generalized  variational  techniques  such  as  Hylleras 
or  configuration-interaction  expansions  using  a  sufficiently 
large  number  of  undetermined  parameters  give  results  for  the 
upper  bound  to  the  energy  to  almost  any  desired  accuracy  for 
small  atomic  systems.   For  example,  the  ground  state  of  helium 
has  been  computed  to  eight  figures  using  a  39  parameter  wave 
function  and  to  ten  figures  using  a   IO78   term  wave  function. 
However,  the  wave  functions  so  obtained  are  not  uniformly 
convergent.   This  is  further  brought  to  light  in  the  fact  that 
the  rigorous  lower  bound  calculated  with  these  wave  functions 
does  not  usually  agree  with  the  upper  bound  in  the  last  few 
figures.   If  the  variational  wave  functions  were  a  close 
approximation  to  the  exact  solution  over  most  of  the  atomic 
configuration  space  one  might  expect  better  agreement  between 
the  two  limits.   Furthermore  the  num.ber  and  type  of  correlation 
terms  which  may  be  Included  in  the  trial  wave  function  are 
usually  chosen  ad-hoc.   Therefore  the  Interpretation  of  rhe 
resulting  variational  wave  function  for  a  two  or  three 
electron  system  does  not  usually  lead  to  any  real  insight 
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into  general  electron-electron  correlation  effects.   Since 
extensions  of  such  expansions  to  atoms  with  many  electrons 
must  become  hopelessly  complex,  some  other  techniques  must 
be  utilized  to  exhibit  general  correlation  effects  and  give 
wave  functions  of  some  reliability  for  larger  atomic  systems. 

In  this  paper  we  develop  the  method  of  asymptotic 
expansions  for  obtaining  atomic  wave  functions.   We  look  for 
asymptotic  large   Z  (atomic  number)  wave  functions  of  the  form 

f  =   e'^'^laQ  +  a^Z"^  +  agZ"^  +•-  +  ).        A 

The  energy  is  also  expanded  as  a  power  series  in  Z   of  the 

form 

p   p    -1         p    _p   p 
E  =  -Z  (Cq   -  Z    C^   +  Z    C2   +••■+)      B 

Insertion  of   A   and  B   into  the  Schrodinger  equation  for  any 
atom  leads  to  a  set  of  lower  order  partial  differential  equations 
for  h,   a„   and   a, •   Although  the  energy  parameters  may  be 

determinable  in  a  self  consistent  manner  we  make  use  of  the 

2         2  5 

fact  that   C„    and   C-,    are  known  exactly  for  atomic  systems. 

2 
Cp   may  be  determined  from  an  experim.ental  analysis  of  the 

binding  energies  of  the  isoelectronlc  series  of  interest. 

In  section  2  of  this  paper  we  apply  our  formalism  to 

separable  atomic  systems  and  show  that  truncated  asymptotic 

expansions  are  exact  solutions  to  the  separable  problem.   In 

section  3  the  partial  differential  equations  for  a  general 

atomic  system  are  solved  recursively  through  the  lowest  orders. 
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Screening  and  correlation  are  exhibited  explicitly  in  the 
resulting  atomic  wave  functions.   The  correlation  effects 
are  exhibited  as  a  sum  of   2   electron  correlation  functions 
which  depend  on  the  principal  quantum  numbers  of  the  electrons 
Involved . 

Asymptotic  solutions  of  form  A  (naturally  with  different 

largeness  parameters )  have  been  obtained  previously  for  a  wide 

6 
class  of  electromagnetic  scattering  problems    and  fluid 

dynamics  problems.   In  those  rare  cases  where  an  exact 

solution  could  be  obtained  by  the  method  of  separation  of 

variables,  a  comparison  between  the  asymptotic  solution  and 

the  exact  one  has  shown:  1)  the  leading  term  or  terms  of  the 

asymptotic  series  represents  the  leading  term  or  terms  in  the 

asymptotic  expansion  of  the  exact  solution.   2)  the  asymptotic 

expansion  including  only  a  few  terms  is  a  very  good  approximation 

to  the  exact  solution  over  most  of  the  domain  when  the  expansion 

parameter  is  of  the  order   2   or  greater.   Prom  (2)  we  expect 

our  asymptotic  expansion  to  be  a  good  approximation  to  the 

exact  solution  even  for  atoms  or  ions  of  low  atomic  number. 

That  this  is  the  case  is  shown  by  the  close  agreement  of  our 

results  of  section  ^  for  the  angular  correlation  in  the  helium 

(Z  =  2)   wave  function  away  from  isolated  singularity  points, 

with  those  obtained  from,  a  many  parameter  variatlonally 

calculated  wave  function.   Prom  (1)  we  expect  our  asymptotic 

series  to  represent  the  leading  terms  in  the  asymptotic  expansion 

in  large   Z   of  the  exact  solution.   This  Is  shown  in  part  at 
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the  end  of  section  4. 

In  section  5  we  calculate  some  first  order  expectation 
values  for  the  helium  ground  state  utilizing  our  asymptotic 
expansions.   Our  helium  wave  function  is  shown  to  obey  the 
virlal  theorem  to  first  order.   Magnetic  susceptibilities 
within  5  percent  of  the  accepted  values  are  obtained  for 
helium  and  singly  ionized  lithium. 

In  section  7  we  determine  by  a  "Matching  Technique" 

the  asymptotic  expansions  for  Helium  which  are  correct  within 

a  small  neighborhood  of  the  singular  points.   The  expansions 

are  shown  to  satisfy  the  "Cusp  Condition"  at  the  origin  exactly 

and  the  one  at  s-.^    (inter-electron  distance)  =0   to  leading 

order.   Since  the  wave  function  expansions  in  these  regions 

£n   Z 
have    Z    terms,  expectation  values  for  operators  other  than 

the  energy  will  include  such  terms  in  their  high  Z   expansions. 

Ordinary  perturbation  theory  expansions  (7a)  (7b)  (7c)  do  not 

include  such  log  terms. 

2.   ASYMPTOTIC  SOLUTIONS  FOR  SEPARABLE  ATOMIC  SYSTEMS 

The  general  Schrodinger  equation  for  a  system  of  a  nucleus 
with  charge   Z   surrounded  by  N  electrons  can  be  written  in 
dlmensionless  form 

N  N  N 

-  I  51  V  2  ^  -  X^f  V'  +  ©5111^   ^  =  E^     (1 
1=1  1=1   1      ^1-1  j>i  ^Ij 

where   s^   is  the  distance  in  Bohr  radii  of  the  1th  electron 
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from  the  nucleus,   E   Is  the  energy  of  the  system  In  large 
Rydbergs  (approximately   27-2   ev)  and  (^  is  set  equal  to 
0   when  the  electron-electron  interaction  is  considered  to  be 
turned  off  and  taken  as  1   for  the  actual  atomic  problem. 
Let  us  consider  the  case  when   (h)  =  0.   The  equation  is 
then  separable.   That  is,   ^  ;ij  (Zf(  s-,  )0(  s^)  •  •  • -(Zil  s^)   where 
^(  s  )   satisfies 

-  I  V^2  0  -  I  0  -  E/  =  0.  (2) 

The  exact  solution  can  then  be  represented  by  product  wave 
functions  of  hydrogenic  character.   If  the  Fauli  principle 
is  invoked  and  the  hydrogenic  states  are  filled  two  at  a  time 
in  the  order  of  increasing  energy,  one  builds  up  the  periodic 
table  with  each  element  having  its  proper  valence  and  the  correct 
ground  state  symmetry  properties  (assuming  L-S   coupling  to  be 
negligible).   Thus  the  independent  electron  model  certainly 
predicts  some  of  the  more  important  atom.ic  properties  correctly 
and  must  in  some  sense  represent  a  first  approximation  to  the 

true  solution.   The  energies  of  such  separable  systems  only 

2 
include  the   Z    term  and  can  therefore  be  written 

2 
E  =  -Z  E   large  Rydbergs.  (3) 

Each  electron  with  a  principal  quantum  number  n   contributes 

o 
an  amount   l/2n   to   E„ .   The  exponential  behavior  of  the  wave 

function  for  an  independent  electron  in  orbital   n   is  given  by 

0^e-^^/^    .  (4) 
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In  the  next  section  when  the  general  problem  Is  solved  we 
will  demand  for   ©  =  0,   C-,^  =  0,   Cg^  =  0,...,  C.  =0 
that  f     go  over  to  a  product  of  separable  wave  functions  each 
of  which  has  an  exponential  behavior  expressed  by  (4).   In 
order  to  make  complete  contact  with  the  known  solution  with 
the  interaction  turned  off  we  ask  that  each  electron  for  the 
given  state  of  the  atom  or  ion  has  the  same  value  of   n  with 
the  interaction  on  or  off. 

In  the  remainder  of  this  section  as  a  simple  Illustration 
of  the  asymptotic  expansion  technique  we  apply  it  to  the 
separable  problem  (2).   Substituting  the  formal  expansion 

00 


j^^^-Zh(s)^   a.(^)Z-J  (5) 

J=0     ^ 


and   (3)   into   (2)   and  taking  the  gradient 

00 


VjZf  =   -ZSZfVh  +   e-^h^I  z"^(Va.) 
and  Laplacian 

00  00 

V^^  =    [Vtf  Yl     a.Z-J+2   _    sY"  {7h.Va,)Z-J+l 

00  00 

-  V^h^a.Z-J+l  +  X^Z-^V^a. 
We  obtain  setting  the  coefficient  of  the   Z^   term  equal  to   0 

-  I  (Vh)2aQ  =   -   EqEq  .  (6) 
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Since   h  =  h(s),   the  equation  can  be  written 


as/ 


'0 


which  admits  the  solution 


h  =  x/2E„'s  =  as 


,1 


Setting  the  coefficient  of  the   Z    terms   =  0   gives 

1      -*■  9  n   "^n      ^n 

-  i(Vh)  a^  +  aQ— ^  +  ^aQ-—  =-  E^a.^  . 

Ss     s       s 


(7) 


The   (Vh)    term  cancels  the   E„   term  and  we  are  left 
with  homogeneous  equations  for  a^^   of  the  form 


^0 


Defining   P  = 


Ss 
1 


-1  -  1 


'0 


=    0. 


(8) 


1^   our  result  for  a^   can  be  written 
a^  0 

a^  =  F(0,E^)s^ 


Setting  the  coefficient  of  the  zeroth  power  of   Z   equal  to 
0,   an  inhomogeneous  equation  for  a,  of  the  form 


aa- 


1     3    2 
P  i^  -  2^  V  ^0 


9) 
a 


1 


results.   In  terms  of  the  variable   r-,   defined  as   r.  =  - —  , 

1  X   a^ 

equation  (9)  becomes 


as 


^V% 


^"'o^o 


(9-1 


which  for  an   S  state  {F{9,0)    =  const.)  takes  the  form 
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^  ^  3(P+l)    _  (9.2) 

ds  ^"^o^ 

Equation   9-2   is    immediately   integrated   and   gives 

r^   =   -   pO+l)/2aQS  (10) 

where  we  have  chosen  the  constant  of  Integration  equal  to 
zero.   Defining  v^  =   ag/aQ  ,   the  equation  for  v^     for  an 
S  state  is 

dr 


2         1     r72 


ds    2aQaQ 


V'a^    .  (11) 


Integrating  we  finally  obtain 

rg  =  0+l)P^0-l)/8aQ^s^  .  (12) 

Thus  to  second  order  in  the  expansion  in  inverse  powers  of 
Z,      ^  is  of  the  form 

!^=.e-^°'0^  s^l  -  z-l^P±^+  ^-2  (P+1)|"(P-1)   ^...^  ).   (13) 

2a„s  8a   s 

For  the   IS   orbital   a  =  1,  ( P=0 )  and  the  series  truncates 
after  the  first  term.  For  the   2S  orbital   a   =  l/2  and 


we  obtain 


0  =  e"^^/^  s(l-2(Zs)"^)  .  (1^) 


which  is  the  exact   2S  wave  function.   In  the  above  manner 
we  can  generate  all  the   S  states  for  one  electron  systems. 
Choosing  F(0,0)   in  a„   as  one  of  the  appropriate  spherical 
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harmonics  and  demanding  truncation,  we  can  generate  the  exact 
wave  functions  for  the  case  when  the  angular  momentum  of  the 
system  is  not  equal  to  zero.  Thus  we  observe  the  well  known 
fact  that  the  exact  solutions  for  the  separable  problem  can 
be  written  as  truncated  asymptotic  series  of  form  (2).   This 

simplicity  arises  because  no  interaction  terms  occur  and  the 

2 

energy  in  the  separable  case  can  be  written   -  Z  E^ =   With 

interaction  terms  present,  the  energy  is  represented  as  an 
infinite  expansion  in  decreasing  powers  of   Z.   This  is  one 
of  the  reasons  no  such  truncation  is  possible  in  the  general 
case . 

3.   ASYMPTOTIC  SOLUTIONS  FOR  GENERAL  ATOMIC  SYSTEMS 

We  assume  an  asymptotic  solution  of  equation  (1)  exists 
and  is  of  the  form 


■Zhi,  s-,,Spj...,s,,  j 


^  =  e     "   "      '^     2^a^{s^,s^,...,r^)    Z   ''.  (15) 

n=0 

The  condition  that  our  solution  go  over  to  the  correct  separable 
form  with  no  interaction  present  indicates  that  there  is  no 
dependence  of  the  function  h   on  angular  coordinates. 
Substituting  equation  (15)  and  the  energy  expansion 

00 


2  \   ,  -,  nW^-w^  2 


E  =  -z^  2_  (-D^'z-^c  ^  (16) 


*  2 

In  the  energy  expansion   C     is  taken  as  necessarily 


positive  for  w  =  0,  I   and   2. 
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into  equation  (1),  we  have  after  cancelling  out  the 

-Zh 
exponential  e    from  both  sides 


oo    N 

-1/2  ZX  [(V'^n^"""'  -  f^V-Vn  +  (vA)a„)Z-"+l 
n=0  1=1  *- 

00     N  00     N 

^   "    J    n=0  1=1  ^i  "        n=0  1=1  J>i  ^IJ 


00 


XI  (-l)"z-^c/5Ia,Z-^  (17) 


=  -  Z^ 


2 

Setting  the  coefficient  of  all   Z    terms  equal  on  both 

sides  of  the  equation  we  obtain  an  eikonal  equation  for  h 
of  the  form 

N 

y  (v.h)2  =  2c/  .  (18) 

1=1 

Since   h   is  not  a  function  of  angles :,  equation  (l8)  can  be 
written 

N 


S  (K 


2       P 

=  2Cq^.  (l8a) 


In  order  that  our  solution  go  over  properly  for  the  case  of 

no  electron-electron  Interaction,   h  must  be  of  the  form 

N 

h  =  x:  fii  (19) 

1=1     1 

where  the  ith  electron  is  taken  to  be  in  a  state  with  principal 
quantum  number  n. .   It  then  follows  that  each  electron  will 

-  1^^  - 


2         2  2 

contribute  an  amount   l/2n    to   C^^  .   Thus   C„   =  E„ 

and  Is  therefore  determined  since   E^^   Is  known.   For  example 

In  the  case  of  a  filled  shell  atom  whose  outer  shell  has  a 

principal  quantum  number  N^, 

N        n-1 

E^  =  2  2II   l/2n^  XI  (  2^+1 )  . 
n=l       £=0 

In   more  conventional  perturbation  theory  techniques  the 
electron  position  coordinates  in  equation  (l)  are  all  multiplied 
by  Z   and  this  results  in  a   l/Z   appearing  as  the  perturbation 
parameter  multiplying  the   l/s .  .   term_o   It  then  follows 
immediately  that  the  leading  term  of  the  energy  expansion  is 

the  unperturbed  energy  obtained  with   l/Z   equal  to   0   which 

2 
is  just  the  separable  energy   -Z  E^^. 

Equating  coefficients  of   Z   terms  on  both  sides  of 

equation  (17)  we  obtain 

N 

y~  [-l/2(^^h)^a^  +  V^heV^aQ  -  l/2(V^^h)aQ 
1=1 

22  ( 20  ) 

-1/s.aQ]  =   -Cq  a^  +  C^  a^  . 

From  equation  18,  it  follows  that  the  coefficients  of  the 
a-,   terms  on  both  sides  of  equation  (20)  are  equal.   The 
resulting  equation  for  a„   is  tnerefore 
N      ^^ 

1  =  1       1       1 
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Equation  21  admits  the  general  solution 


where 


N    -CO.  A.s. 
a,  =J^TT^s,    ^e^-  (22) 


o).  =  1  -  n. 


^1  =  i^iCi^/N 

and   J_   Is  an  arbitrary  function  of  angles  and  of  the 
differences  formed  by  taking  the  radial  coordinates  of  any 
two  electrons,  multiplying  each  by  Its  principal  quantum  number 
and  then  subtracting  one  from  the  other.   Choosing   J^^   as  a 
product  of  the  appropriate  spherical  harmonic  functions  of  the 
angular  coordinates  of  those  electrons  not  In   S  states,  our 
solution  goes  over  correctly  for  the  separable  case.   From  the 
above  form  of  the   A.'s,   it  follows  that  in  any  given  ion  or 
atom  the  shielding  factor  for  any  electron  is  approximately 
proportional  to  its  average  distance  from  the  nucleus  with  the 
interaction  turned  off. 

Equating  coefficients  of  the  zeroth  power  terms,  we 
obtain  an  inhomogeneous  equation  for  a.,   of  the  form 

ii     ^  +    ih     -  1)^  a,)  -  C^a. 


H,  BiT^n,  -  "-'I.    ^l'    -    ^1  ^1 


2„     ^  2, 


1=1   -1  --1     -1      -1 

N  N 

-  ®  a^  ZIZ:  ^^.  +  121  V.%  -  C^'a^  .   (23) 
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Making  the  substitution  a-,  =  ^i^n      ^^'^   noting  that  both  functions 
a^  and   a,   satisfy  the  same  homogeneous  equation^,  we  find 
equation  23  reduces  to 


<    -.    n.  ds  . 

1=1   1    1 


N 


N 


-~-®^l 


+ 


1=1  j>l   ij 


kl 


0  1=1 


^1% 


24) 


The  solution  to  equation  24  can  be  written 
N  N 

-J 


^1    -    ©   \Ylc..   +|lIrK 
1=1    j>i 


/^V.    2^ 
1     a 


1=1 


0    ds. 


a 


0 
N      p    2 


1=1 


N      ^1^1 


(25) 


where      C.  .      is    the   solution   of 
ij 


1     Sc     1   1     Sc 

1 

n.    ds  .         n .    ds  . 

s 

1         1           J         J 

IJ 


(26) 


subject  to  the  boundary  condition  that  as  either  s.   or 
s  .   goes  to  infinity  or  zero  with  the  other  held  f Ixed^, 

C.  .   goes  to  zero.   This  condition  follows  since   C.  .  represents 

ij  ij 

the  direct  correlation  term  in  our  first  order  solution.   As  the 
electrons  go  infinitely  far  apart  from  one  another,  the  electron- 
electron  Interaction  contribution  to  the  wave  function  goes  to 
zero.   The  solution  of  equation  26  is  shown  by  standard 

o 
techniques    to  be  of  the  form 


ij 


■n  .n. 


's  .  .Q.  .  +  (  n  .-D.  .n.  )s  .  +  (n.  -D.  .n  .  )s  . 
1,1  1,1      3      1.1  1   1      1   ij  J   J 


A. 


27) 


ij 
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where 


0.  .   is  the  angle  subtended  at  the  nucleus  between 
electrons   i   and   J,   D. .  =  cos  0. ., 


2      2 
Q-,- -;  =  \/n.   +  n.  -  2D.. n.n.,   and 

A.  .  =  (Q.  .+n  .-D.  .n.  )(s  .-(n  ./n.  )s  .)   for  n  .s  .  <  n.s  .  . 

It  is  obvious  from  2?  that   C.  .   depends  only  on  the  ratio 
s . /s  .   and  the  angle  between  the  electrons. 

It  should  be  noted  that  in  the  analysis  of  this  section 
we  have  tacitly  assumed  that  no  non-uniformities  occur  in 
performing  the  two  limiting  procedures  (w)   — >  0,   Z  — ^   <». 

4.   ASYMPTOTIC  GROUND  STATE  SOLUTIONS  FOR   2   ELECTRON  ATOMS 

The  ground  state  of  a  two  electron  atomic  system  consists 

2 
of  a   Is    configuration.   Since  both  electrons  are  in  orbitals 

with  a  principal  quantum  number  of   1^   we  obtain 
CO.,  =  cOp  =  0 

o 

A-,  =  Ap  =  A  =  C-,  /2  , 


and  equation  22  for  a   gives 


A(s-,+Sp) 
aQ  =  e   ^   ^   .  (28) 


Operating  on  a^  with  the  Laplacian  we  obtai 


n 
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V.%    2A 

— ^ =  +  A 

Performing  the  Integration  in  equation  25  we  are  led  to 
the  result 

^1  "  C-^^{1,1)    +   A  log  s^s^  +  1/2(,A^-C2^)(s^+S2)       (29) 


where   C.  .{n.,n.)      refers  to  the  correlation  function  when 
electron   1   is  in  an  orbital  with  principal  quantum  number 
n.   and  electron   j   is  in  an  orbital  with  principal  quantum 
number  n..   Equation  27  for  n.  =  n.  =  1   gives 

-,       /s,  ^  +  d(s-,+Sp)  \ 

C,p(l.,x)  =  -  i^  log  ( -^^^ ^--^  (30) 

^  V(l^d)|s^-S2l   / 

with 

d  =  (l-D-,^)/^  and  D-,^  =  cos  e.^- 

The  isoelectronic  energy  expansion  is  well  known  for  the 

case  of  2  electron  atomic  systems.-^   In  terms  of  large  rydbergs 

we  observe  that   C.   =  .625  and   Cp  =  .I5765.   Using  these 

values  we  obtain  A  =  -3125  and  the  coefficient   p   of  the 

S-,  +  Sp   term  is   -  .  03  •   The  first  order  asymptotic  wave 

function  for  a   2   electron  atomic  system  is  then 

-(Z-A)(s.+Sp) 
^^   =  e        ^     [1  +  Z  ^(0^2(1,1)  r  A  log(s^S2}+P(s^+S2)  )]" 

(31) 
Since  the  value   .3125  for  our  shielding  parameter  A   is 
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exactly  the  value  found  by  performing  a  variational  calculation 
on  a  simple  exponential  product  wave  function,  the  leading 
term  of  f^      is  already  well  known.   On  graph  1,  C-^^''^'-^'^ 
is  plotted  as  a  function  of  x  (x=S2/s-|_)   for  G^g  =  ^    '      "^^^^ 
is  the  case  when  the  two  electrons  are  both  on  the  same  side 
of  the  nucleus  and  lie  on  a  straight  line  which  passes  through 
it.   We  observe  that  C.^     goes  to  zero  as  the  electrons  go 
infinitely  far  apart.   As  the  electrons  come  near  one  another, 
C-,  p  becomes  more  and  more  negative  leading  to  a  smaller 
value  for  ^.   and  therefore  to  a  decreased  probability  of 
finding  the  electrons  close  together.   These  are  the  correlation 
effects  one  expects  physically. 

In  point  of  fact,  C-.^     diverges  logarithmically  at 
S-,  =  Sp.   For  the  case  of  Q-.^  =   0°   the  divergence  is  stronger 
and  a  first  order  pole  occurs.   This  type  of  singularity  at  a 
point  about  which  the  asymptotic  expansion  changes  character 
is  common  to  many  asymptotic  solutions.   For  exam.ple  the 
WKB  asymptotic  solution  of  the  sim.ple  harmonic  oscillator 
diverges  in  the  neighborhood  of  the  classical  turning  point. 
The  asymptotic  solution  also  changes  character  about  this 
point.   Often  one  can  obtain  by  a  stretching  technique  an 
accurate  non-divergent  asymptotic  solution  in  the  transition 
region  where  the  original  asymptotic  solutions  diverged.   For 
the  harmonic  oscillator  such  a  stretched  solution  is  the  Airy 
function.   Simply  stated,  such  singularity  phenomena  occur 
because  the  postulated  form  of  the  asymptotic  expansion  is  not 
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adequate  to  represent  the  asymptotic  expansion  of  the  true 
solution  In  this  region  due  to  problems  of  non-uniform 
convergence.   This  will  be  discussed  further  In  sections  6 
and  7  when  we  obtain  the  asymptotic  expansion  In  the  n.b.h. 
of  singular  points.   In  another  sense  one  may  sometimes 
think  of  the  occurence  of  such  singularities  as  due  to  our 
conversion  of  a  problem  Involving  a  second  order  partial 
differential  equation  to  one  where  a  set  of  first  order 
equations  are  solved  recursively.   Somewhere  in  the  process 
a  boundary  condition  may  be  lost.   In  trying  to  recover  the 
correct  boundary  value  the  asymptotic  solution  overshoots  the 
mark.   Thus  we  do  not  expect  our  asymptotic  expansion  to  be 
valid  in  some  neighborhood  of  the  points   s-,  =  Sp.   However 
from  the  plots  of  our  correlation  function  as  a  function  of 
X  as  shown  In  figure   2  for  various  angles,  we  expect 
radial  correlation  to  be  a  general  property  of  the  helium- 
ground  state  system.   That  is,  one  should  observe  a  reduced 
probability  of  finding   2   electrons  having  the  same  value 
of  the  radial  coordinate.   Strong  radial  correlation  effects 
in  the  helium  wave  function  have  been  discussed  elsewhere. 
We  also  observe  from  figure  2  that  the  correlation  effect 
for  any  value  of   x  decreases  with  decreasing  D-,  ^  and 
therefore  with  increasing  s-,p  as  expected. 

In  figure  3,  we  compare   ^-,   for  helium  and  doubly 
ionized  beryllium  for  the  case   s.,  =  -5^  0-,  ^  =  0  .   It  is 
physically  unreasonable  that  the  true  helium  wave  function 
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becomes  negative  In  the  neighborhood  of   s-j^^  =  '^-   ^^  note 
from  the  figure  that  this  phenomenon  will  occur  much  further 
from  the  point   s.  =  s^   in  the  case  of  helium  than  for  doubly 
ionized  beryllium.   Therefore  it  is  immediately  obvious  that 
our  solutions  become  applicable  closer  to  s^  =  ^2  ^®   ^ 
becomes  larger.   As   Z   goes  toward  infinity  our  solution 
^,   will  become  applicable  over  almost  all  configuration  space. 

It  should  also  be  pointed  out  that  IS  orbltals  give 
rise  to  a  logarithmic  singularity  when  either  electron  goes 
into  the  nucleus.   The  region  in  the  neighborhood  of  the  nucleus 
in  which  this  effect  is  large  is  small,  even  for  Z  =  2.   It  is 
of  the  order  of   .3   bohr  radii  for  this  case.   A  non-divergent 
asymptotic  solution  which  is  correct  in  the  neighborhood  of 
the  origin  can  be  found  by  a  stretching  technique.    As  our 
electrons  go  an  infinite  distance  away  from  the  nucleus,  the 
usual  non-convergence  problems  arise  in  the  wave  function. 
Since  this  domain  will  not  be  of  interest  in  computing 
expectation  values  we  do  not  trouble  ourselves  about  it  here. 
To  avoid  this  problem  one  has  to  use  a  density  matrix  formulation 
of  the  problem. 

In  figure  4,  we  plot  V-i   as  a  function  of   Sp  for  the 
case   S-,  =  2,0,  p  =  0  .   In  the  same  figure  we  plot  the  leading 
term  of  i/-.    which  is  just  a  product  of  screened  exponentials. 
We  refer  to  this  term  as  ip^.      Note  that  f^      becomes  less 
than  ^Q   in  the  neighborhood  of   s-,p  =  0  and  larger  than 
f        as  the  distance  between  electrons  increases.   Because  of 
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the  compensating  effects  in  the  different  regions,  when  a 
correct  asymptotic  description  of  the  wave  function  In  the 
Immediate  neighborhood  of   s,„  =  0   and  the  origin  Is  used 
In  conjunction  with  our  asymptotic  solution  which  Is  valid 
everywhere  else,  we  expect  the  correct  normalization  constant 
multiplying  ^-,   will  be  only  slightly  changed  from  the 
normalization  constant  for  ^_ .   Assuming  tp^      to  have  the 
same  normalization  factor  as  f^,      we  have  calculated  a 
correlation  correction  factor  which  Is  just  the  ratio  of  the 
two  wave  functions.   In  figure  6,  we  compare  our  results 
with  those  of  Gilbert   who  obtained  this  correction  factor 
utilizing  a  many  parameter  varlatlonally  determined  Hylleras 
type  wave  function.   In  these  plots  we  are  holding  the  radial 
coordinates   s.,   and   Sp  fixed  while  changing  the  distance 
between  the  electrons  by  varying  the  angle  subtended  at  the 
nucleus  by  the  two  electrons    (for  0-,  ^  =  0  ,   s-,^  =  s.-Spj 
for  ©ip  =  l8o°,  s-|2  =  s-|+Sp.)   We  expect  physically  the 
greatest  correlation  to  occur  for  0.  p  smallest.   This  is 
borne  out  by  the  curves.   In  the  case  of  the   s-,  =2,  Sp  =  3 
case  we  run  Into  the   S-,   approximately  equal  to   Sp 
distortion  phenomena.   Our  correlation  factor  Is  too  low  for 
small  angles  as  expected.   However  outside  of  this  region  our 
results  agree  quite  well  with  Gilbert's  and  If  we  used  the 
correct  normalization  constant  for  f-,      we  expect  that  the 
agreement  would  be  Improved. 

For  large   s.,   and   Sp   It  Is  known  that  the  exact  solution 
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for  helium  has  an  exponential  behavior  given  by 

Tl,  ^  Q-  /^   (S1+S2)  (32) 

where   E   is  the  energy  eigenvalue.   Putting  in  the  energy 
expansion  for  E   in   terms  of  the   C's  and  expanding  the 
exponential  to  first  order  in  powers  of   l/Z  we  obtain  for 

Z  large 

-{Z-A)(s-,+Sp)       .1       o        o 
f  ^   e  ^"^[l  +  Z  ^{A  -C2  )(s^+S2)/2]. 

Comparing  our  asymptotic  expansion,  equation  2$,    with  this, 
we  observe  that  it  gives  the  correct  high  Z   asymptotic 
expansion  of  the  exponential  part  of  the  exact  wave  function 
for  the  above  case. 

5-   EXPECTATION  VALUES  FOR   2   ELECTRON  ATOMIC  SYSTEMS 

In  the  previous  section  we  investigated  the  properties 
of  the  asymptotic  helium  wave  function  and  found  certain 
unphyslcal  logarithmic  divergences  which  occur  at   s.  =  0, 
and   s.  =  s..   However  these  inaccuracies  in  our  wave  function 
have  little  effect  on  the  expectation  values  of  those  operators 
which  are  relatively  Insensitive  to  the  values  of  the  wave 

function  in  the  neighborhood  of  these  points.   Such  operators 

2  2    ? 

as   s^2  '      ^12'   l/2(s-j_  +S2  )>  1/2(3^+32),   and   1/2(1/3^+1/32) 

fall  into  the  above  class.   Expectation  value  computations  to 
first  order  in   l/z  were  performed  for  these  operators  and 

our  results  and  a  comparison  with  the  "exact"  results  of 

ii 
Pekerls   for  helium  and  singly  ionized  lithium  are  given  in 


2k 


the  tables  below.   All  values  are  glA^en  in  atomic  units 


Table  1.   Helium  Expectation  Values* 


Operator 

Expectation  Values 

Pekeris 

Asymptotic 

l/2(s-^^+S2^) 
l/2(s^+S2) 

l/2(l/s-^+l/s2) 

2 

^12 

^12 

1.193 

.9295 
1.68832 
2.516 
1.^22 

1.186 

.9455 
1.68750 
2.600 

1.547 

+ 


Table  2.    Lithium  Expectation  Values 


Operator 

Expectation  Values 

Pekeris 

Asymptotic 

l/2(s^^+S2^) 

.4463 

.4624 

l/2(s^+S2) 

.5728 

.5868 

1/2(1/3^+1/32) 

2.6879 

2.6875 

2 

3^2 

.9271 

.9925 

3^2 

.8623 

^9307 

Our  calculations  necessitated  one  numerical  double 
integration  over  a  finite  domain  for  each  expectation  value ■ 
These  integrations  were  performed  on  the  7094  computer  usin| 


Pekeris '  data  is  originally  given  to   10   figures 
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a  trapezoidal  rule  numerical  integration  scheme  .   Our 
numerical  integration  results  are  generally  good  to  4 
figures.   For  the  calculation  of  the  nuclear  magnetic  shield- 
ing constant,   l/2{ l/s^+l/s^) .   the  results  were  extrapolated 
to  six  significant  figures.   The  value  of   1.68750   for  this 
expectation  value  shows  that  Vj   satisfies  the  virlal  theorem 
to  first  order  in  l/Z."*"^  Our  result  for  the  diamagnetic 
susceptibility,   l/2( s^^+Sg^) .  of  helium  is  within  l/2  per 
cent  of  the  exact  one. 

In  performing  our  computations  we  form  the  product 
ip     Otp-.      where  0   is  the  operator  under  consideration  and 
calculate 

0  = ^ (33) 

to  first  order  in  l/Z.   When  performing  integrations  not 
involving  terms  dependent  on  s  p   or  D-,  ^   it  is  convenient 
to  choose 

00  00 

[  dV  =  v^tt)  J  J  ds-,dSpS-,  Sp  . 
0  0 

For  integrations  involving  correlation  terms  we  take 

CO   3        11  "1  "5 

TdV  =  27r^/dsf  du  uf  dt(s^-t^)  (3^) 

0   0,0 

where   s  =  s-,  +  s^j   t  =  s.,  -  Sp  and   u  =  s-,p.   For  operators 

of  the  form  0.  =  l/2(  s-j^'^+S2"^ ) ,   the   A(  logs-j_+logS2)   term  is 

evaluated  utilizing 


26 


I 

0 


,    -kx  j-i 
dx  e    x'^log  X 


,1.'[psl(,1+l)-log  k] 


k 


J+1 


(35) 


where 


dq 


psi(J+l)  =  T  dy  e'^y-^log  y  =  ^  log  P(  q  )| 
0 

14 


'j+1 


and  has  been  evaluated  elsewhere.     Making  a  change  of 
variables  in  the  integral  containing  the   C-,  p   term 


u  ->  a  =  u/s 

,,  ,  2  ,2.1/2/,,  2  ,2vl/2  ^  a/2 
t^y=  (l-(u  -t  )    /((s  -t  )    a)) 

the   s   integration  can  be  performed  separately  and  what 
remains  is  a  double  integration  over  a  finite  domain  of 
the  form 


J_  -L 

CJ  =  2jda  a(l-a^)Jdy  f  .(y,a)(log  W(y,a)  -  log  y)   (37) 
0  0     ^ 


with 


W 


l_(l_y2)a)(2-y^) 
(l+(l-y2)a) 


1/2 


and 


f  (y,a)  =  (0  /s 


3,/((l_(l./,V)5/2(2-y2)l/2, 


The   log  y   term  in  equation  37   can  be  evaulated  numerically 
by  rewriting  it 


fdy  f  (y,a)log  y  =  -  f  .(0,a)  + 
"O     ^  ^       0 


dy(f  .(y,a)-f  .(O,a))log  y.   (38) 
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We   obtain  finally 


0.  =  (iilil 


j         2(2Z)'^ 


fl  +  |[(j  +  2[psl(3+j)-psl(3)]) 
.  1  5i   CO  -  i6±izl4     cjI  +  0(1/Z 


2. 


(2+j)!2 


(39) 


Similarly  expectation  values  of  operators   s-,^   are 
found  to  be  of  the  form 


J  _  (  J+5  )  •'   /J:. 


12    (2Z)'^'25  (  j+3   3(j+5) 


1  + 


I  (a(j  +  4[psi(6+j)-psi(3)]  -  2  log  4 


2  I 


1 


-'(+  ^^  CO 
2  2-^ 


CS12J 


(40) 


3+j    3(5+j 


_1 1 

3+j   3(5+j) 


where   CS12J  has  the  same  form  as   CJ   except  that   0  ./s"^ 

J 

Is  replaced  by  a'^   In  the  Integrand  and  we  have  taken 

1       .  1 
I.  =  f  dPp^"^-^'  /  da(l-a2p2)log(l-a2p2) 
J    0        0 

which  can  be  evaluated  analytically  in  terms  of  standard 

quadratures . 

6.   Technique  of  "Matched  Asymptotic  Expansions" 

The  asymptotic  expansion  for  the  ground  state  of  helium 
has  been  found  to  diverge  at  certain  singular  points.   We  shall 
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now  discuss  the  technique  of  "matched  asymptotic  expansions" 
(stretching)  for  obtaining  an  accurate  non-divergent  asymptotic 
solution  in  the  neighborhood  of  any  one  of  the  singular  points. 
Let  us  consider  the  exact  solution  ^(  Z.;x-,  ^x^^x:^ )   of  our 
many  body  problem  with  Hamiltonlan  H  and  a  two  term  formal 
asymptotic  expansion  f-.      obtained  by  the  methods  of  section  5  ■■ 
Suppose  we  want  to  determine  the  asymptotic  expansion  of  ip 
in  the  neighborhood  of  the  point   x-,  =  0   where  f-.    blows  up 
unphysically .   Now  although  it  has  not  been  generally  proven, 
we  take 

9^(Z  large)  — >  ^^ 

where  the  wiggle  over  the  f      indicates  "the  expansion  of" 
and  the  arguments  Inside  the  parenthesis  indicate   the 
limiting  procedures  for  obtaining  the  expansion.   If  more 
than   1   limiting  process  is  involved,  the  one  furthest  to 
the  right  is  the  first  to  be  performed.   Now  when  we  evaluate 
Tp-,      as   X.,  — ^  0   we  are  performing  the  limit  process 

^(x^— >0,Z  large)  -^  f-^{x-^-^G)   ->  divergence. 

However  near  the  point   x-^  =  0  for   Z   large  out  finite  we 
really  want  to  evaluate 

V^(Z  large, x-j^->0) 
where 

^(Z  large, x^->0)  ^  ^(x^-^,Z  large)  ^  Tp-^ixj-^O)    . 
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Let  us  consider  a  transformation  of  variables  where  we 
rewrite  our  operator  H   In  terms  of  a  new  stretched  variable 
i-..      That  Is 

X  — >  ^-1  =  Z^   X,   where   j  >  0   ( j   as  yet  otherwise 

unspecified ) 


K[Z;x-^,x^,x^)    -^  H{Z;i-^,x^,x^)  . 

Such  a  change  of  variables  guarantees  that  when  considering 
the   Z   dependence  of  terms  in  H: 

A.  Terms  in  H  which  go  as   x-,  (k>0)   will  be  small  when 
expressed  in  H.   On  the  other  hand   k  <  0  terms  will  be 
large.   [Note  that  for  x.,   small  we  expect  terms  ^  x-^    (k>0) 
in  H   to  be  small . ] 

B.  Terms   -^   in   H  will  be  of  magnitude  Z"^  -^   in  H. 


in  the 


This  implies  that  if  we  look  for  a  solution  of  the  form 
^  R5  0  exp  f  (  Z,  ^-1  ,Xo^x^  ) ,   the  leading  term  of  -^ 

L         ci.        J  OX-, 

neighborhood  of   x.,  ;5^  0  will  be  of  magnitude  Z"^   greater 

than  the  leading  terms  of  -v^   or  -v^   in  this  region  when 

ox  2      ox^ 

each  is  reexpressed  as  a  function  of   ^-,  ,  Xp  and   x^ .  [Note 
again  that  for  x-,   small  we  expect  our  boundary  layer  solution 
to  vary  more  rapidly  with  x    (have  a  larger  partial  derivative) 
than  in  the  case  of  the  asymptotic  solution  away  from  the 
singular  point.   The  proper  stretching  accomplishes  this]. 
Thus  we  see  that  stretching  is  a  way  of  ordering  the 
terms  in  H   in  the   n.b.h.  of   x-,^0.   When  considering 
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H,   this  ordering  is  expressed  simply  by  the  dependence  of  the 
various  terms  on  a  single  parameter  Z.   We  can  again  look  for 
solutions  of  the  standard  asymptotic  form  (15)-   The  resulting 
system  of  partial  differential  equations  to  be  solved  is 
naturally  simpler  than  the  original  problem.   However  because 
of  our  ordering  the  equations  are  only  applicable  near  x,  =  0. 
Let  us  call  the  formal  two  term  asymptotic  solution  to  the 
stretched  problem   -  f-^{Z;^-^,x^,x^).      We  expect  that  if   j   has 
been  chosen  correctly  that  "ip-.      is  an  asymptotic  expansion  of 
^  which  is  valid  for  Z   large,   s-,   small.   Since   ^-,  =  Z'^x  , 
the  solution 

^■^(^^  large) 

corresponds  to  f{x-.    small,  Z  large)   and  therefore  we  expect 

^^(l-j^  large)  — >^-j^(x^  small). 

This  is  called  the  "matching  condition".   Only  the  correct 
value  of   j   will  allow  this  condition  to  be  satisfied  even 
though  certain  constants  and/or  functions  in  f-.      are 
determined  from  it.   Once  it  is  satisfied,  the  solution   ^-, 
for   I-,   small  obviously  corresponds  to  f{Z   large, x.,  small) 
and  is  therefore  the  desired  asjmiptotic  expansion  in  the 
neighborhood  of   x.  ^0.   Another  important  consideration  in 
determining  the  requisite  value  of   j   is  that  we  expect  at 
least  some  of  the  operators  present  in  the  P.D.E.   for   r-, 
of  the  original  asymptotic  expansion  will  also  be  present  i 


n 
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the  P.D.E.   for  m-,  .   In  addition  we  also  expect  some  second 
order  partial  derivative  terms  to  appear  in  the  equation  for 
m.,  .   For  the  case  of  ip^      becoming  divergent  at  x-,  =  0, 
Xg  =  0   it  may  be  necessary  to  perform  different  stretchings 
on  variables   x-,   and  Xp  but  essentially  the  same  consider- 
ations go  through  for  understanding  the  "matching  principle." 
From  physical  considerations  one  can  usually  ascertain  whether 
f     varies  more  rapidly  (its  partial  derivative  is  greater) 
with  respect  to  one  variable  than  the  other  in  the  neighborhood 
of  the  singular  point.   It  is  then  the  former  variable  which 
receives  the  greater  stretching. 

In  the  forgoing  we  have  introduced  the  idea  of  stretching 
and  have  indicated  how  the  matching  principle  is  fulfilled. 

We  now  state  more  explicitly  the  form  of  the -matching 

(15) 
principle     utilized  in  this  paper  which  is 

"The  two  term  expansion  of  tj/..       (when  it  is  expressed 

in  terms  of  stretched  variables)  in  powers  of   /z  =  the  two 

term  expansion  of  f      (when  it  is  expressed  in  terms  of 

ordinary  variables)  in  powers  of   /z." 

For  the  case  that  ^-,   is  singular  at  x-,  =  0,   defining 

^^(Z;^^,X2,x^)  =  ^^(Z;— ,X2,x^)   , 

flE   "  j^i(  EXPANDED  FOR   Z   LARGE,  KEEPING  2  LARGEST  ORDERS), 

'y 

^(Z;x^,X2,x^)  =  ~^(Z;Z"^x^,X2,x^)  , 

SlE  "  ^i( EXPANDED  FOR   Z   LARGE,  KEEPING  2  HIGHEST  ORDERS) 
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and 


the  "matching  principle"  asserts 

7.       SOLUTIONS  IN  THE  NEIGHBORHOOD  OP  SINGULAR  POINTS 

Before  we  begin  our  investigation;,  we  would  like  to 
summarize  the  notation*  used  in  this  section. 

f-^,    The  helium  ground  state  asymptotic  wave  function  obtained 
in  section  4  is  written  in  terms  of  ordinary  variables  as 

Tjy-.,      When   ^-,   above  is  rewritten  in  terms  of  the  stretched 
variables  associated  with  the  singular  point  under  investigation 
it  is  denoted  by 

/      -Zhr,  ,  1    -I 

ti  =  '     --^^  ^  Z^^ 
where   r-,   is  written  in  terms  of   d,,  .   In  general  a  wiggle 
under  a  symbol  means  it  is  being  reexpressed  in  terms  of 
unnatural  variables . 

il/^„,      ^-i  „   is  the  expansion  of  1//-.       in  powers  of    /Z,    keeping 
only  the  two  leading  terms. 

r-,  „,   r-,_   is  the  expansion  of   r.,   in  powers  of    /Z , 


f-^. 


keeping  only  the  leading  term. 

^    The  asymptotic  solution   of  the  stretched  problem  in  terms 
of  stretched  variables  is  denoted  by 

^  =  e    [1  +  2" J 
*    In  region   lb,   these  definitions  are  modified  slightly. 
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where   h'  =  h  for  regions   I  and  they  are  not  equivalent 

in  regions   II. 

't}/,    When  ^  is  written  in  terms  of  unstretched  variables,  we 

denote  it  by 

-7-    -Zh'r,  ^  ^l^ 
f  =   e      ^   [1    +  Tr-\  • 


ii/_,  ili^   is  the  expansion  of  ^   in  powers  of   /Z,   keeping 
only  the  two  leading  terms. 

m-i  ,   m-,    is  the  expansion  of  m-,   in  powers  of   /Z, 
^Ig   ^Ig  ^1 

keeping  only  the  leading  term. 

Our  helium  ground  state  asymptotic  solution  ^-,   diverges 
in  the  following  regions: 
la.    (t=0,  Ut^O) 
lb.    {u=0,  ST^O) 
Ila.   {s-^=0,  SgT^O) 
lib.   (s^=0,  3^=0) 

To  find  the  correct  asymptotic  solutions  for  regions  I  we 
rewrite  the  2  electron  Schrodinger  equation  in  coordinates 
s,  t,  u  as 

,(32_t2)(^  +  ^  +  ^)  +  2s(u2-t2)  tl^ 
Ss    Su    St  SsSu 

+  2t(s2-u2)^  +  4  su$^  +  2(s2-t2)|i  -  4  ut|^         (41) 

aust  ^  ^^  ^ 

+  4Zsu^  -  (s^-t^)^  +  Eu(s^-t^)^  =  0 
with  E  given  as  before  by   {l6).   Regions   la   and   lb 
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are  considered  removed  from  the  origin.   Thus  we  do  not  have 
to  consider  any  stretching  in   s   and  expect  that  the 
exponential  dependence  on   s(  shielding  effect)   will  "be  the 
same  for  the   la   and   lb   solutions  as  for  the  original 
asymptotic  solution.   Therefore  we  look  in  regions   la  and 
lb   for  stretched  solutions  of  the  form 

^  ^  e-(Z-A)s^-^  ^  ^  ^  ^  ^._^]/  (42) 

In  appendix   I  we  have  substituted  the  above  into  equation  4l 
and  have  kept  terms  through  order   /z   considering  %i  -  ili{^,t,\i) 
Certain  terms  automatically  cancel  one  another  in  equation  Al = 
These  terms  are  enclosed  in  boxes  and  need  not  be  considered 
further.   Now  consider  the  general  double  stretching 

^  ^  Z^t 
7  =  Z"^u. 

Under  each  term  of  equation  Al   we  express  the  order  of  that 
term  in  Z   when  considered  as  a  function  of   4,  j,      and  Z 
in  terms  of  the   j,  k   of  the  transformation o   For  the  case 
j  =  0,  k  =  0   (no   stretching)   the  largest  terms  are  of 
order  Z°   and  equating  them  equal  to   0   we  obtain 

us  [A  -2^^]  -  ut  [A  -2-^] 


Sm 


-    2s[u  -t  ]-^  +  4asu  -  s   +  t   -  C^  [us  -ut  ] 


0. 


It  is  a  simple  matter  to  show  that  this  equation  is  satisfied 


*      f      here  is  taken  as  the  general  expansion,  but  only  term.s 
through  m-,   will  I 
"matching  theorem". 


through  m-,   will  be  of  importance  in  satisfying  the 
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by  the   r,   (equation  29)  of  our  original  asymptotic  expansion. 
We  will  Indicate  this  as  expansion  A.   The  terms  in  equation 
Al   entering  the  leading  equation  for  expansion  A  (equation  42) 

are  indicated  by  writing  an  A   above  the  term.   Since  there 

2 
is  one  term,   -  s    in  Al   which  is  Independent  of   j   and   k, 

it  is  necessary  that  other  terms  in  Al   which  depend  on   j 

and   k  are  also  of  order  Z°   if  we  are  to  be  able  to  equate 

the  coefficient  of  this  power  of   Z   equal  to   0.   For  the 

region   la   no  stretching  is  required  on  the  s-.^     variable 

and   j  =  0,  k  =  /2  are  the  proper  values  of   j   and   k. 

For  the  case   j  =  1,  k  =  /2  we  obtain  the  expansion  in 

region  lb.   For  both  cases  the  P.D-E.  to  be  satisfied  by 

m-,   is  found  by  setting  the  coefficient  of  the   Z    terms 

equal  to  zero.   The  terms  entering  the  equation  for  m-,   in 

region   la   have  a   la  written  above  them  in  eq .   Al . 

Similarly  those  entering  in  region   lb   have  a   lb   written 

above  them.   Note  that  of  the  four  differential  operator 

terms  present  in  A,   two  are  present  in   la   and  one  is 

present  in   lb.   If  for  example  we  had  chosen   j  =  1,   k  =  1 

as  the  proper  stretching  in   lb,   none  of  the  differential 

operator  terms  present  in  A   would  have  entered  the  equation 

for  m-,  .   Needless  to  say,  pursuance  of  the   j  =  1,  k  =  1 

solution  for  m-,   would  have  indicated  that  no  match  with  the 

asymptotic  expansion  A  was  possible. 

We  now  write  down  the  P.D.E.   for  m.,   in  region   la   as 
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B  m       Sm  ^m  u  op 

1/2 
where   ^  =  Z  ■^  t.   Now  we  have   r-,   of  expansion   A   which 

can  be  written 

2   2 
r,(s,t,u)  -  [Alog(^-^)+  i[A2-Cp^]s-  ^log(^  +  ^^   )] 
^  i|      2     2     2d     (i+d)lt| 


lA 


45) 


2 
2  ^2    "^ 


where   d  =  [— ^ — n]         ■>   This  is  Immediately  seen  by  referring 
s  -t 

to  equation  29  and  rewriting  it  in  terms  of  the  variables 
s,   t  and   u.   Following  our  prescription  at  the  end  of 
section  6 

2      2 

Zs' 


2        d 
r^    =    [Alog  ^  [1  -  ^^]    +   V2[A^-C2^]: 


^^     (1+dJ  Vz^/^l 

with  .2      1/2 

1   -   ^  A,  2 


£  1  j 

d^(Z;|,s,u)  =  d(Z;  Jy^  ,s,u)  =  -[ g ^ 


Zs^ 


Expanding  for  large   Z 


.2 

s~    u" 


l^s  ^  u[l  +  |2(^  -  ^)  +  0(^)] 


we  obtain 
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'SlE 


J  2 

Alog^  +   V2[?.^-c/]s    +  |j[log  ^j^ 

La 


(46) 


L 


-    (log   2u-log(l   +   Vs))]l  +   0(i^    ) 


In  order  to  satisfy  the  matching  principle  m, „   obtained 
from  m-i   of  equation  44   must  be  equal  to   r,  „  given 
by  equation  46.   Choosing 

m^  =Alogs^  +  ^/2[A^-C2^]s  -  s/2u[  log2u-log(l  +  Vs  )  +  M-^  ] 

(47) 
and  substituting  Into  equation  46  we  find   M,   satisfies  the 
heat  flow  equation 


-  2  — - 


ae 


2 


(48) 


To  complete  the  match  we  expect  for  large   ^ 


M  =  -  log  ^  +  ^  a.(s)^  '^      with  a  I  .1  =  0.       (49) 


Substituting  49  Into  48  we  obtain 


^^S 


J+2 


j-2)( j-i)aj_2^ 


aa. 

2  -^  2  J 


0 


(50) 


Thus   a^  =  constant.   Since  ^  must  be  even  In   t  and 
therefore  in   ^   we  next  set  the  coefficients  of   |~^   In 
(50)  equal  to   0   obtaining 


2 


+  1  =  0 


(51) 
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which  has  the  solution 

B.^   =  ^/2  s.  (52) 

For  higher  values  of   j   we  find  recursively 

aj  =Lt2Hiii  ps  aj_,   .  (53) 

O 

Defining  w  =  s/|   we  observe  that  for   ^   large,   w  small, 
we  can  write 

M^  =  -  log  ^  +  a^  +  ^  +  ^  w^  +• • •+   .     (5^) 
Therefore  we  look  for  a  general  solution  of  (^8)  of  the  form 

M^  =  -  log  ^  +  a^  +  0(w)  (55) 


with 


!Z^(w)  -^  w/2  +  5/4  w^  +•  ■  •+  (55a) 


for  w   large.   Substituting  we  obtain 
,2, 


d' 


:^+  1  r3  _  1  1  M.  _  ^_  (56) 


^  2  ^  2  ^w     2J  dw      h„2 
dw  w  Hw 

which  has  the  general  solution 

1  1   ^     -^ 

,  r.   -^/2  "2^    if.   "^^2  "  2^r^  .  e^^'       /.^^ 

0  =  C  J  dww  /  e     -  I]-  J  dww    e     J  dw  ■  — —,  '    ^  ^7 ) 


w'  ^2 


substituting  V  =  p—  we  obtain 


1  ^   V^/2  ^  "^  V  /2      J  u'^/2 


V^/2 
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Now  the  homogeneous  term  would  admit  odd  powers  of   ^   which 
is  unphysical.   Therefore  we  set   C-,  =  0.   Now  we  write  our 


solution 


(t)  =   1  f^  e-\(V)  (59a) 


with 


(V)  =/  ^  e^  .  (59b) 


The  lower  limit  on  integral  5913  has  been  chosen  so  that 
no  contributions  from  partial  integrations  will  result  from 
it.   For  V  large,   w  small  we  partially  integrate  to  higher 
negative  powers  of  V  giving 

g(V)  ;i;e^[-^  +  ^/2  4-^+  2^   !_+...  +  ]       (6o) 
y-^/2        V^^  2         2     V''^2 

which  results  in 

^,,,,1^  li5_^_i£.^  i_        ,61, 

n=l      (2    )n      V 

which  gives  as  its  leading  terms 

(t)(V)  ;^  i^  +  ^-p  +-•+  =  I  +  I-  w^  +•'•+   .         (62) 

i6v 

which  is  the  previously  obtained  expansion  (55a). 
We  can  now  write 

m^    =Alog  s^   +   ^/2[A^-C2^]s    -   ^/2u[  log2u-loga+Vs  ) -log|+aQ 

+  (t)(w)].  (63) 
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We  rewrite   w  =  / //^   =  ^/^^^  ,   Thus  for   Z   large   w   is 
small  and   (t)(w)   is  given  by  (62).   Therefore 


m^  ^Alogs   +  ^/2[A^-C2^]s  -  ^/2u[  log2u-log(l+Vs  ) -log^+aQ  ] 

-  \^-n   +  O(^).       (64) 
^  4ut^     Z^ 

Truncating  64  according  to  the  prescription  at  the  end  of 

section  6  to  obtain  m,  _  we  note  that  for  m.,  ^   to  match 

A/IE  ~1E 

r-i  „     we  must   choose 

Hq   =   ^/2   log   Z.  (65) 

Now  that  we  have  determined  M-,   completely,  to  obtain 

the  expansion  near   ^  =  0  we  partially  integrate  59t>  to 

increasing  powers  of   V   obtaining 

1  2 

g(V)  =   e^[-  -f-  +  ^V  /2(i    2  ^  _  2   v2  +...  +  )j     (gg) 

y-i/2  ^  ^-^ 

which  leads  to 


(1)(V)  :i^  ^  log  V-V  +  ^  ? ^^ —  .   (67) 

2  j  =  2   [l-3--(2j-l)](j) 


In  terms  of  our  variables   ^   and   s 

2       2       4 

2 
(1)  »loge  -  |log2s  -  11^  (68) 

o 

Substituting  into  (55)  for   ^  /s   small  we  obtain 

2 
M^  R>  -  I  log2s  +  \   logZ   -  I  ^  /s         (69) 
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and  our  solution  m,   no  longer  shows  the  singularity  for 

t  small  which  appeared  in  the  original  asymptotic  expansion 

r-,  •   Thus  the  correct  solution  to  be  used  in  region  la  is 

2 
s   ,  lr.2  ^  2- 


1^  =  A  log  -^  +  ■^['h    -Cg  ]s 


+  |^[-  I  logZ  -  log2u  +  logd  +  ^)+  |log2s 

2        2  2 
+  If^   +  0(f^  )  ].        (70) 


Now  let  us  consider  region  lb.   We  choose  the  stretching 

V2 

7  =  Zu  . 
Therefore 

2  .2 


3J1   =   e-^^-^)^[l   +  |[Alog  |-(1    -  ^)    +l/2[-K^-C^^]s  +  C^^]] 

(71) 


"  ■  Zs 


with 


^^^^      ^<^a^°«^,i.d,,lV,V,| 


and 


(71a) 


For   large      Z 


2.1'^     1/      J^   ^    2zL    2         .2J-' 
Z    /2    s  s  ^ 


^v/ 


leading   to 


'v. 


212-^"    ''{'^Tfr+    V2[^+    1]  (72) 
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Truncating  (keeping  3  terms  this  time)  we  find 

llSZ   ^  11  W  nfJ. 


Now  from  the  appendix  we  see  that  the  equation  for  m.   in 
region  lb  takes  the  forni 

S  m^      £  S  m.,      A  2   -,  \  Sm-i    -, 

Using  a  power  series  construct  as  In  case  la  we  are  led  to 
an  m-j   of  the  form 

m^   =  aQ(s)    +   /Xiw)  (75) 

.2 
where      w  =  /s      as   before.      Substituting   (75)    Into    ( 7^ ) 

we   find    that     3C(w)      satisfies 

4w  ^+   2[w  +   UX    -   1   =   0  (76) 

which  has  the  solution 

w 

Y=C-f—  e      +  _J__  e       J   4—  e  /2   ^w.    (77) 


W 


^/2  4w^/2         -^   w^y2 


In  order  that  our  solution  be  even  in   t  we  choose   C  =  0. 
For  w   large  we  Integrate  by  parts  to  obtain 

w 
For  Z   large,   w   large,  comparison  with  (73)  Indicates  that 
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one  should  take 

2 

Eq  =  A  log  V^  +  V2[A^-C2^]s  +  V2       (79) 

1   "^/2 
in  order  to  match  everything  except  the   /z  '    term. 

Now  for   t  small,   w  small  one  can  Integrate  (77)  to 

higher  powers  of   w   to  obtain 

;^(w)  =  V2  yZ  (-i)'' ^"^   •  (80) 

n=0  3  •  5  •  •  •  (  2n+l ) 

1/2 
We  now  note  that  if  we  had  Included  an  m-,   /Z  ^    term  in 

V2 
expansion  ^2,      m,  /   would  have  satisfied  the  homogeneous 

1/2 

part  of   equ.  (7^)  and  any  function  m.,  /  ( s ,  | )   is  an 

1/2 


acceptable  solution.   We  shall  now  show  that  m.,  /   can  not 

1/2 

enter  our  asymptotic  solution  in  the  neighborhood  of  ^^0. 
Let  us  first  write  down  our  asymptotic  solution  in  region 
lb  as 

^^e-(^-^)^[l+(l/2ju(l  -i(^)  +0(^)') 

(81) 
"^1/  ^^)  ^^  2 

+  T +  ¥'^    l°g  ^  A  +  V2(A^-Cp^)s  +  V2)] 

Z  /2       ^  "^ 

From  part  B  of  section  6  we  expect  that  in  region  lb 

^U,7,s)  ^   0(Z^/2)  ||(^^^^3)_  (82) 


One  observes  that 
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'dm 


M^e-(2-A)Sf_.^2^^^  J^ 


Bt 


and 


12s 


]/2. 


'2    at 


^    12s  ^tt; 


1!=^ 


Z-A)s    1 


am 


i^    at 


^. 


83a) 


(83b) 


Now  assume   m  /  (s  ? )   Is  expressible  as  a  power  series  about 


V 


2 


^  =  0.   For  any  power   j 


am 


^ 


i^^zJ/V-1  =  j|^'-^z'/2 


.0 


a 


n 


84) 


in  addition 


This  would  contribute  a  term  of  order  Z    m 

to  the  term  already  present  of  order  Z  ^    and  condition  (82) 

would  not  then  be  fulfilled.   Thus  in  the  vicinity  of  ^-^0, 

m-,     must  vanish  identically  and  therefore  will  not  contribute 

to  (8l).   Thus  finally  for  the  solution  in  region  lb  we  write 


f   «  e-^^-^^^[l  +  1/2  u(l  -  i(^)  +  0(^)^) 


2 


(8la) 


+  h7^    log  ^  A  +  V2(A^-Cp^)s  +  V2)] 


2 


.ti6) 


Now  according  to  a  theorem  proved  by  Kato'^^S   the 
exact  wave  function  satisfies  the  cusp  condition  at  u  =  0 
of  the  form 


au 


=  V2 


(82) 


u=0,t=0 


45  - 


Performing  these  operations  on  our  solution  8la  we  find 
that 


^u 


=  [1  - 

u=0 
t=0 


|(Alog  ^  A  +   s  +  ^/2)] 


(83) 


Therefore  in  the  limit  that   Z  — ^  °o   this  cusp  condition  is 
identically  satisfied. 

In  order  to  obtain  our  expansion  for  one  or  both  of 
the  electrons  near  the  origin,  it  is  convenient  to  rewrite 
the  Schrodlnger  2  electron  Equation  in  terms  of  variables 

ds-j     s-|  ds-,   dSp    So  dSg      du 

I       2    2 ,2>  ^2,      ,   2    2,  2>  ^2, 
h      b-il/         ^^1  ~^2  +u  )  d  ^      [s^    -Sn^    +u  )  d  ^ 
+  —  —    -^   

du 


S-,u     Ss-,Su 


SgU 


SspSu 


+  2(E  +  —  +  ^ 

S-,     Sp    u 


)^  =  0  . 


(84) 


For  the  case  that  electron  1  is  near  the  nucleus,  let  us 
assume  an  integral  stretching  of  the  form  i-,    =   Zs-,    is 


correct.   Noting  that  in  terms  of   u,  s, ,   and 
in  r   is  expressed  as 

u^-(Sp-s,)   ^^2 

4s-|^S2 


d 


(85) 


-  46  - 


we  can  write 


^T  =  e 


wlth 


and 


(z-A)s    -i    Ae  /    r  i 


e    e 


+  V2(A2-c/)(^-f  s^)  +  C^2^^ 


/vl2 


2d. 


log 


U  +  d  ( 75—  +  R  ^ ) 


(1 


d,  = 


u^  -  (s. 


i^)(s2 


/Z 


'!/„) 


2  11 


Z    2 


/Z 
/2 


(86) 


(86a) 


For  larse   Z 


72 


'w  1'*"' 


'l'^2 


72 


{86b) 


which  leads  to  an  expansion  for   Z   large  of  the  form 


42 '^ 


^  tu 


'2  J 


+  O(^) 

z 


(87) 


Therefore  truncating  to  two  orders 


<s/ J.- 


e    [1  +  |-(-A^nZ  +  A^^ 


(88) 


Ain^^  +  A  insg  +  ^ /2[-h^-C^]s^)] 


This  suggests  that  we  look  for  a  solution  of  the  stretched 
equation  of  the  form 


-  ^7  - 


+  ^  nip  +••'+).        (89) 
Z 

Making  this  ansatz  and  substituting  into  the  stretched  P.D.E. 
B^l      |,  S^      SSp     Sp  ^Sp      ^u     u  Su 


>2,  2 


li         _  2^   2    d^   +  2   —    '^ 


(90) 


+  ^^ z  ^ ^-^^  +  z 

u    S^-,  Su       |-,u  S^i^u      I-,   d|-,Su 

u    Z  uSp    Sp  SspBu 

2 
+  2( -Z^( Cq^-Z~^C^^+Z"^C2^+" ' •+)    +  ~  + -O^  =  0 

we  obtain  from  setting  the  coefficient  of   Z    terms  equal 


to   0 


a^Ji)       p      •■    am,(i) 

^  +  [—  -  2]— i—  =  -  c/  .  (91) 


A  solution  of  this  equation  is 


2$ 

1   .^  „^1 


24,     ^  J       e 


m|"'  =  \[i-^+!lni^-   -~-  -   I  f   — ^"2—  +  C]         (92) 


'1         -00       ^1 

where   C   is  a  constant  to  be  determined.. 

Prom  setting  the  coefficient  of  the   Z      terms  equal 
to  0   we  obtain 

-  48  - 


;^  (2) 
-  2[— ^ ^  +  l[c/-A^]] 


+  1  +  2[^  -  1]    2 


(95) 

+  -l-[u^-s/][— ^  +  ^]    -  §  =  0. 


4-,u         du    SuS^ 


1 


Choosing  to  set  the  first  parenthesis  itself  equal  to  0 
we  obtain 

For  Z   large  we  can  rewrite 

Zs,       2Zs. 


m5^^-;A[Zs.+inZs,-  -^ —   f   ds   ^ p-+C] 

-^^  ^  ^      2Zs^    2Z   -00    ^   s^"" 


92a) 


and  performing  the  usual  truncation  procedure 

i|g^-A[Zs^+inZs^+0(|)  +  C]  .  (92b) 


Zi 


Considering   (94)  and  (92b)  in  89  we  observe  that  the  matching 
condition  with  (88)  is  made  if  we  choose 

C  =  -  ^nZ.  (95) 

Thus  the  matched  stretched  solution  is 


49  - 


(,Z-AjSp        ^-j  -|  _  p  p 

^  =   e  e      ^[1   +  |[A^ns2+    /2[A^-C2    Jsg 


+   Mi-^+^n^^- 


2?- 


-00  ^ 


(96) 


To   obtain   the   expansion  for     s-,      small,    we    Integrate  by  parts 


to   higher  powers    of      4      to   obtain 
f        e 


e      ^[1   +  ~[--K£nZ+l\£ns^+  ^/2(A   -C^    )s2 


A(-i-+   0(|,^))j] 
3  ^ 


and  therefore 


_    -(Z-A)Sp  -Zs^      .  (Zs,  )  ^ 

f   =   e  ^e    -^[1  -  f( ^--+0(Zs^)^) 

+  |[-AinZ+Ains2+   ^/2( A^-C^^ )s2 ] 


97) 


for  the  case  of   s-,   small  (In  region  Ila )  .   Thus  the  correct 
asymptotic  expansion  has  no  divergences  at   s.,  =  0.   The 
solution  shows  the  expected  physical  characteristic  that  for 
one  electron  very  near  to  the  origin  the  shielding  effect  of 
the  outer  electron  upon  It  goes  to   0. 

The  cusp  condition  at   s-,  =  0   satisfied   exactly  by  the 
correct  wave  function  Is 


Ss 


1 


(98) 


s,  =  0 
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Now  for  our  asymptotic  solution  (97) 


Ss-, 


7^+   -^^-^^"2  -2^1  -2ZAS, 

2^  +  ^         ^     [ I+o(s2)] 

3  ^ 


(99) 


which  in  the  limit   s-,  =  0  becomes 


as. 


ZV'(S2=0) 


100) 


s^  =  0 


Our  asymptotic  solution  identically  satisfies  this  cusp 
condition. 

In  section  la  we  perform  additional  integral  stretching 
on  Sp  and  u.  We  define  i-.  =  Zs-,  ,  ^p  =  ZSp  and  7  =  Zu. 
One  can  readily  show  that  for  this  case 


-( ?  +^  ) 


with 


and 


+  MJini^+£ni^)    +   C^2(^^,^2'^^^^ 


'^12^1^  ^2'  "^   ~    "2^ 


2  /^   /.  ^2^   ^/2 
7  -(?2"^1^ 


1  ,.  /z:!!ii:kL.) 

2^      V(i+d)U2-^il/ 


41 


'1^2 


101) 


Transforming  equation  84  to  new  coordinates  we  obtain 
the  usual  "perturbation  theory"  equation  in  the  form 


51 


(102) 


i!i_  +  2_  23i^  ^  h^  +  ^^]L.  +    2  ^ 
S^/    ?3_  S?l    a^2^    ^2  ^^2      ^^^ 


+  2((-Cq^+Z~^C^2-Z"^C2^+--'  +  )  +  |-  +  l-)^'  "I  f  ^  ==^  0 
We  look  for  a  solution  of  102  of  the  form 


+   ^  m^{i^,^^,n)+-'-+^   .       (103) 

Substituting  (103)  into  (102)  one  finds  that  m^^  and 

(  2) 
m-j^    satisfy  equations  of  the  form  (91)  in  their  respective 

variables.   Therefore  from  the  analysis  for  section  Ila   it 

follows  that  m-j     and  m-j     will  be  matched  with  the  first 

three  terms  of  the  coefficient  of   /Z   in  equation  (101). 

We  find  the  remaining  equation  satisfied  by  M-,   to  be 

^  +  2[^  -  1]  _J.  +  1  +  2[1-  -  1]-  1 


ai^  ^1  a^        s^2  ^2    ^^ 


2 


+  2  —^  +  7  -^  -  [^ — ^— ^  +  — — ]  — ^ 

by  '^  ^7        ^^7         4^7       ^7 
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+  [_i 2 ]  — l-H— ]  — —=      f  .        (104) 


Since  we  have  already  taken  proper  care  of  the  functions 
depending  only  on  either  variable   ^-,  or  ^^      near  the  origin 
with  our  choice  of   m-j     and   m-j     we  expect  that  near  the 


orisin 


M^  =  C7  +  0( 7£)  (105) 


where   C   Is  a  constant  to  be  determined  and  £     represents 
one  of  the  fundamental  smallness  distances  7,    i-.      or  ^^■. 
Substitution  of  (105)  into  ( 104 )  gives   C  =   /2.   Although 
we  have  not  shown  the  matching  condition  fulfilled  for  M-, 
with   C.  p   we  have  determined  the  leading  term  of   M-,   near 
the  origin.   Thus  for  region  lib  keeping  only  the  leading 
terms  of  m-j    and   m     we  have 


^^  e"''le"''2[l-2Alf  +(l/2Ju  -  ^is^^^s/)    .   0(|1)] 


(io6) 

assuming  no  other  constants  are  brought  in  by  matching  M^ . 
Note  that   106   identically  satisfies  the  cusp  condition  (98) 
for   S-,   or   Sp   equal  to   0.   The  cusp  condition  at   u  =  0 
is  exactly  satisfied  in  the  limit  that   Z  — >  <»  •   Thus  in 
the  neighborhood  of  the  origin 


^  « 1  -  Z(s 


^+s^)    +   ^/2u  -  2A^  +  0{f)  (106a) 
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(17) 
in  the  limit  of  infinite   Z   our  f     gives  correct 

behavior.   This  is  to  be  compared  with  perburbation  theory 

where  the  coefficient  of  the   u   term  near  the  origin  does 

not  go  over  to  its  exact  theoretical  value  in  the  infinite 

/  -1  O  \ 

Z   case     .   Since  ordinary  perturbation  theory  assumes  at 
the  outset  no   log  Z   terms  are  present  in  the  wave  function 
expansion  whereas  we  have  shown  that  in  certain  spaclal 
regions  they  must  be  included,  it  is  not  surprising  that 
the  perturbation  theory  wave  functions  and  expectation 
values  calculated  from  them  will  show  certain  fallings. 
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SPATIAL  CORRELATION    S,  =  .5,   ^,2=0° 


Figu  re  3 
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ANGULAR  CORRELATION-HELIUM 
_L    vs  ai2 

Gilbert's  Correlation  Fa ci^r. 

Case  of  S|  =  I.O,  $2=  3.0  (RMP  6/63)^ 
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